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The neutron-rich isotopes with Z<20, in particular those with neutron numbers around N=28, 
have been at the focus of a lot experimental and theoretical scrutiny during the leist few years. 
Shell-model calculations using the effective interaction SDPF-NR, were able to predict or to explain 
most of the properties featured by these nuclei. Prominent among them is the disappearance of the 
N=28 shell closure for Z<16. We have incorporated into SDPF-NR some modifications, either on 
purely theoretical grounds or guided by new experimental information. The proposed interaction 
SDPF-U offers enhanced reliability with respect to the earlier version. 
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I. INTRODUCTION 

In this article wc describe the reshaping of the effec- 
tive interaction SDPF-NR, frequently used for Ohu shell 
model calculations in the sd-pf valence space. This in- 
teraction originates in the work of Retamosa et al. [l|, Q 
and its final form was fixed in More details about it 
and some applications can be found in and in many 

of the experimental articles cited below. From the very 
beginning, the interest in the region was fostered by the 
experimental indications of a weakening (or even a com- 
plete disappearance) of the N=28 closure [1,0, Hi and by 
the predictions of several mean-field calculations in the 
same sense [H, [ll|, [l^l . How can a shell closure vanish? 
Its persistence depends on the balance of two opposite 
tendencies. On the one hand, magic numbers are associ- 
ated with energy gaps in the spherical mean field and this 
means that to promote particles above a closed shell costs 
energy. On the other hand, this energy can be partly re- 
covered, because closed shells have no correlation energy, 
while open-shell configurations of neutrons and protons 
have a lot. Several examples of this phenomenon exist in 
stable nuclei in the form of coexisting spherical, deformed 
and superdeformed states in a very narrow energy range. 
A fully worked out case study in doubly magic **'Ca can 
be found in [l3|. At the neutron-rich edge, the struc- 
ture of the spherical mean field may be at variance with 
the usual one at the stability line. The reason is that, 
at the stability line, the T=0 channel of the nucleon- 
nucleon interaction has a stronger weight relative to the 
T=l channel than it has when the neutron excess is very 
large. When the gaps get reduced, open-shell configura- 
tions, usually two-neutron excitations across the neutron 
closure, take advantage of the availability of open-shell 
protons to build highly correlated states that are more 
bound than the closed-shell configuration. Then the shell 
closure is said to have vanished. '^^Mg is a classical ref- 
erence case for this kind of physical behavior. 

The interaction SDPF-NR predicts a clear breaking of 



the N=28 magicity already at S. The closed-shell con- 
figuration (0/7/2)*^ represents only 24% in "48, 28% in 
^^Si, and just 3% in *°Mg H. Its description of all the 
sulfur isotopes is excellent |14l | and its prediction of an 
isomer 0"'' state jus t above the first 2+ has been verified 
experimentally 15*. In addition, according to the cal- 
culations, the ground state of ^'^S should be | instead 

of the expected | , suggesting that this nucleus is de- 
formed. This has been verified experimentally as well 
[H, [l3l. Many other properties are also in agreement 
with the experimental data [l^, [l^ . 



The main reason that has prompted us to improve 
the existing effective interaction has been the vigorous 
renewal of interest in the region due to the access to 
many new neutron rich species, in particular silicon iso- 
topes, and indeed, because of the conflicting experimen- 
tal evidence on the nature of ''^^Si cor ning from MSU and 
GANIL 0, [U, m, m, m, m, [la 113, p. some 
of these experimental data will be analyzed in the final 
part of the article. As more and more experimental re- 
sults are accumulated, the need for a modified version of 
the SDPF-NR interaction becomes more compelling. 



As advanced in the title, our valence space comprises 
the sd shell for the valence protons, covering from Z=8 
to Z=20 and the sd and pf shell for the neutrons, that 
is from N=8 to N=40. Notice however that below N=20 
these calculations reduce to the pure sc?-shell ones, and 
beyond N=Z=20 to pure p/-shell ones. The lightest nu- 
cleus in the valence space, that plays the role of inert 
core, is therefore ^^O. The heaviest one is ^°Ca. The ef- 
fective interaction is composed of three blocks; sd-sd, sd- 
pf (proton- neutron), and pf-pf (neutron-neutron). All 
the shell-model calculations are unrestricted in the full sd 
shell for the protons and the full pf shell for the neutrons. 
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II. THE REFERENCE STATES 

In order to improve the monopole part of the SDPF- 
NR interaction, we rely on the new experimental infor- 
mation coming from nuclei that have one particle or one 
hole on the top of neutron and proton closures, as well 
as on the values of the quasi-particlc gaps across these 
closures. Indeed, the correlations arc fully taken into 
account. This means that in order to find the right 
monopole changes several iterations are usually needed. 
The results for the reference states, using the new inter- 
action SDPF-U, are compared in Table U with the exper- 
imental data (taken from the NNDC database [1^ when 
no explicit reference is given below) and with the predic- 
tions of its precursor interaction, SDPF-NR. The neutron 
gap in ^®Ca is 4.74 MeV, compared with the experimental 
value, 4.78 MeV. 

The nucleus '^^Si plays a pivotal role in fixing the evolu- 
tion of the effective single-particle energies in this space. 
The reason is twofold; it is very single-particle-like, be- 
cause of the doubly-magic character of '^^Si, and its neu- 
tron single-particle states define the location of the pf 
orbits at the neutron-rich side. That is why the changes 
made in the initial interaction of Rctamosa et al. [l| to 
comply with the experimental finding of the low excita- 
tion energy of the | state of "^'^Si Q had such important 
consequences in some cases. The modifications leading to 
the interaction SDPF-NR left the calculated controls of 
the i and | states of ^^Si at their initial positions, 
i.e. similar to the ones they have in ^"'^Ca. Besides that, 
in order to place the | state of '^^Si at its experimen- 
tal excitation energy, only the Qd^/2-^P?,/2 centroid was 
modified. As a consequence, SDPF-NR produces a ^ 

ground state in ''^K, instead of the |^ favored by the 
experiments [29l | , and docs not reproduce the assignment 
of J=0^ to the ground state of proposed in [30] ■ 

Two more pieces of experimental information are con- 
sidered as well; the recent study of the single-particle 
states in '''^Ar [23| that has shown a certain reduction 
of the spin-orbit splitting of the p orbits compared with 
the situation in *^Ca, and the spectrum of ^''S [sH. Fi- 
nally, the cross-shell monopoles were also constrained by 
requiring exact agreement with the pf-shell centroids in 
''^Ca extracted from the ■*°Ca(d,p) reaction (3^ . 

With all these boundary conditions, the monopole part 
of the interaction would be uniquely defined, provided we 
knew the full ■^^Si single-particle spectrum, which we do 
not. Among the possible choices, the one that fits most 
naturally with the available experimental indications [33j 
is that of a compressed spectrum with excitation ener- 
gies 1 MeV, 2 MeV and 4 MeV for the |", \~ , and |" 
states. This time, we produce the compression by chang- 
ing equally the 0(i5/2-r3, 0^3/2-13 and lsi/2-r3 centroids 
(rp is a shorthand notation for all the orbits in the major 
oscillator shell of energy hio(j)+ |) except the one with 
the largest j = p + 1/2). In order to obtain reasonable 
binding-energy differences, we have adjusted the global 



TABLE I: Predictions of the SDPF-U interaction for the ex- 
citation energies of the reference states, compared with the 
experimental data and with the predictions of the SDPF-NR 
interaction (in MeV). 
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French-Bansal monopole parameter a=^(3V -|-V ) of the 
sd-pf interaction to the one-neutron separation energy of 
^^Ca (8.35 MeV) and the total centroid of the neutron- 
neutron pf-pf interaction to the ^^Ca - ^°Ca mass differ- 
ence. This translates into a modification of all the T=:0 
and T=l sd-pf centroids of -0.045 MeV and a modifica- 
tion of aU the T=l pf-pf centroids of -0.250 MeV. If the 
experimental data on separation energies would demand 
it, extra terms quadratic in the number of sd valence 
protons. n^=(Z-8), and p/ valence neutrons. n^=(N-20), 
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FIG. 1: (Color online) Effective neutron single-particle ener- 
gies at N=28 from Z=8 to Z=20 . 



FIG. 2: (Color online) Effective proton single-particle energies 
in the calcium isotopes from N=20 to N=40 . 



could be added to the interaction without changing its 
spectroscopic properties 

The spherical mean field that is produced by the new 
interaction is better visualized by the evolution of the ef- 
fective single-particle energies (ESPE) with varying pro- 
ton and neutron numbers. These are presented in Figs.[T] 
and m where we plot the p/-shell neutron ESPE's at 
N=28 and Z=8-20, and the sd-shell proton ESPE's at 
Z=20 and N=20-40, respectively. In Fig. [TJ we can fol- 
low the reduction of the N=28 neutron gap as we ap- 
proach the neutron-rich edge. This reduction, combined 
with the degeneracy of the 0d3/2 and lsi/2 proton orbits 
at N=28, shown in Fig. [21 is at the origin of the vanish- 
ing of the N=28 magicity for Z=16. The reason is that 
when the Od3/2 and lsj^/2 arc degenerate, they form a 
pseudo-SU3 doublet that enhances the quadrupolc cor- 
relations of the configurations with open-shell neutrons. 
Notice also that beyond N=28 the orbits cross again re- 
covering their standard ordering. The monopole part of 
the SDPF-U interaction (without the two global modi- 
fications discussed above) is collected in Table fVl at the 
end of the paper, and compared with the monopole part 
of the SDPF-NR interaction. 



III. PAIRING AND THE SILICON CHAIN 

The excitation energies of the first 2+ states of 
the N=22 isotoncs (two neutrons in the p/-shell) are 
predicted too high by the SDPF-NR interaction. The 
deviations increase as protons are removed from *^Ca 
tiU 36Si. The 2+ states of ^^Si and ''°Si are also too 
high by about 400 keV (this fact as been also discussed 
recently in [1^). However, the same interaction gives 
the right excitation energy for the very neutron-rich 
sulfur isotopes. What is the origin of this discrepancy? 
In our opinion, it is a subtle point of effective interaction 



theory. The p/-shell effective interaction contained in 
SDPF-NR is renormalized in order to compensate for 
the absence of 2p-2h excitations from the core. This 
renormalization. that is predominantly pairing-like, and 
that seems to have the right size for Z>14, depends on 
the value of the "^^Ca gap. However, when we deal with 
the silicon isotopes, the relevant gap for a perturbative 
estimation of the pairing renormalization is the splitting 
of the 0(^5/2 and O/7/2 proton orbits in ^''Si, which may 
double the former one. In addition, the off-diagonal 
pairing matrix element {Ody^{J = 0)|I^|0/^^2(^ = 0)) 
is larger than {Odl^^^J ^ 0)\V\Of^/^{J = 0)). Using 
the values of the gaps and these matrix elements, we 
have estimated the excess of pairing renormalization 
for Z<14 to be w250 keV. Therefore, we have re- 
moved from SDPF-NR a schematic pairing Hamiltonian 
in the p/-shell with G=-G.09 MeV, which, for the 
(0/2/2(J=0,T=l)|I/|0/2/2(J=G, T=l)) matrix element 

represents roughly 250 keV. With this change, the 2+ 
excitation energies of the silicon isotopes agree nicely 
with the experiment. Therefore, we have been led, and 
for good reasons, to have two interactions one for Z>14 
and another for Z<14. The files with the two versions 
of the SDPF-U interaction will be available with the 
interaction package attached to the open version of the 
code ANTOINE, that can be downloaded from the URL 



http : / /sbgat 194 . in2p3 . f r/~theory/axLtoine/inain . html 



IV. THE Z=14 GAP; DEFORMATION IN "^SI 

As we mentioned in the introduction, there has been 
recently an interesting debate on the nature of ^^Si, fol- 
lowing the claim of its magic character in Ref. [131 . In a 
subsequent paper by the same authors (2^ , the claim was 
essentially suspended as conflicting evidence appeared. A 
gamma ray of very low energy observed at GANIL was 
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TABLE II: Properties of Si for different effective interac- 
tions. Energies in MeV, B(E2)'s in e^ fm* and Q's in e fm^. 





SDPF-NR 
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SDPF-U 


Z=14 gap 


6.08 


6.22 


5.36 
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N=28 gap 


3.50 


3.20 


2.75 
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-52 
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including tfie 


pairing correction discussed in 


tfie text. 



attributed to the decay of the first excited 2+ state. Such 
a low-lying 2+ state is totally incompatible with a magic 
structure in this mass region [3, H^]. The interaction 
SDPF-NR, aheady indicated the vanishing of the N=28 
shell closure and the development of oblate collectivity 
in "^^Si. Its predictions for ^^Si are gathered in the first 
entry of Table [TTl Even if the excitation energy of the 
2+ state is not low enough, the quadrupole transitions 
are large and consistent with the predicted (oblate) spec- 
troscopic quadrupole moment. Besides, the occupancy, 
6.16, of the Ofr neutron orbit is very low, corresponding 
to a closed-shell percentage of 28%. This evidence was 
somehow masked by the too strong p/-shell pairing that 
we have discussed in the previous section, as can be seen 
clearly in the column of Table M labeled SDPF-NR(") 
(the SDPF-NR interaction including the pairing modifi- 
cation mentioned before) . Even if the closed-shell proba- 
bility increases, the quadrupole collectivity remains and 
the 2"*" excitation energy is now much lower. The other 
difference between the SDPF-NR interaction and the new 
SDPF-U interaction is related to the evolution of the 
Z=14 proton gap from N=20 to N=28. In fact, with 
the interaction SDPF-NR, the centroid of the (proton) 
Ods spectroscopic strength in ''^K is too high in energy 
compared with the experimental value extracted in [s^ . 
The effect of reducing this gap is to lower the 2+ excita- 
tion energy and to increase the quadrupole collectivity. 
Eventually, the SDPF-U interaction makes ''^Si a well- 
deformed oblate rotor with /3~0.4. 

In Table Hill we have gathered the basic spectroscopic 
predictions of the new effective interaction for some 
very neutron rich silicon isotopes and compared them 
with the available experimental data. The quadrupole 
properties are computed with the effective charges 
qi,=0.35, qi/=1.35 that give the best reproduction of the 
quadrupole properties of the sd-shell nuclei [35j. The 
agreement with the available data is very satisfactory. 
The two most neutron-rich isotopes show very distinct 
and peculiar behaviors. In ^^Si, the yrast sequence 
0+, 2+, 4+ does not follow the rotational J(J-l-l) law. 
However, its quadrupole properties (B(E2)'s and spec- 



TABLE III: Basic spectroscopy of the silicon isotopes. Ener- 
gies are in MeV, B(E2)'s in e^ fm'* and Q's in e fm^. 
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troscopic quadrupole moments) arc consistent with the 
existence of an intrinsically-deformed oblate structure. 
In addition, the calculation produces a low-lying excited 
0+ state at 1.08 MeV. This state, the second-excited 2+ 
state, and the second-excited 4+ state are connected by 
E2 transitions of the of the same strength as those of the 
yrast band, but this time, the spectroscopic quadrupole 
moments have similar absolute values but opposite sign. 
This is a blatant case of coexistence of two deformed in- 
trinsic states, one oblate and another prolate. In "'"Si, 
the B(E2) values of the yrast band are smaller than the 
corresponding ones in *^Si, while still being substantial. 
The calculation produces an excited 0+ state that lies 
higher (1.82 MeV) and a very low-lying second-excited 
2+ state at 1.25 MeV. This 2+ is connected by a strong 
E2 transition to the first 3+ state at 1.8 MeV, forming a 
sort of incipient 7 band. 

V. THE VERY HEAVY MAGNESIUMS 

The spectroscopy of the magnesium isotopes has un- 
veiled many cases of unexpected physical behavior. The 
N=Z member of the chain, ^''Mg, is the paradigm of a 
well-deformed nucleus that can be described in the labo- 
ratory frame by means of Elliott's SU3 model. At N=20, 
•^^Mg is not a semi-magic nucleus but, on the contrary, a 
well-deformed prolate rotor. It provides one of the first 
examples of the vanishing of a magic closure by the ac- 
tion of deformed intruders [3^, [13, [11] . As can be seen in 
tab le IIV[ t he scarce experimental information on '^^Mg 
[l§,IMl4l|,[43l and ^^Mg [H] also point to well-deformed 
cases. We produce a '^^Mg less collective than the exper- 
iment, surely due to the absence of neutron excitations 
from the sc?-shcll to the p/-shell in our valence space. In 
■^^Mg, the Ohu! calculation already overshoots the low- 
ering of the 2+ state and, therefore, it docs not seem 
necessary to resort to a large intruder mixing to explain 
it. When more neutrons are added, ■^^Mg, and, promi- 
nently ^°Mg, are extremely deformed prolate rotors with 
very low 2+ states. ^°Mg is at the edge of the neutron 
drip line, and, according to our calculations, it has nearly 
two neutrons in the Ipa orbit, a situation that may favor 

the formation of a neutron halo; ^°Mg could then be the 
first well deformed nucleus adorned with a neutron halo. 
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TABLE IV: Basic spectroscopy of the Magnesium isotopes. 
Energies are in MeV, B(E2)'s in e^ fm" and Q's in e 
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FIG. 3: (Color online) Energy splitting (in MeV) between 
the lowest i and | states in the neutron-rich isotopes of 
potassium, chlorine, and phosphorus . 

In addition, the calculations suggest that in '^''Mg, the 
prolate yrast band may coexist with an oblate one based 
in an excited O"*" state at 2.07 MeV, in parallel with the 
situation found in '*^Si. -^^Mg is, according to our descrip- 
tion, a triaxial rotor in view of the presence of a 7 band 
based upon the second 2+ state at 1.07 MeV excitation 
energy. This band comprises the first 3"*" and the second 
4+ states. The energies and the ratios of the transitions 
between the two bands are consistent with 7=0.2. Sim- 
ilar triaxial structures have been recently shown in the 
heaviest Argon isotopes . 



VI. THE EVOLUTION OF THE 
SPLITTING IN THE K, CL AND P ISOTOPES 

In Fig. [3l we present the evolution of the energy split- 
ting between the lowest ^ and | states in the neutron- 
rich isotopes of potassium, chlorine and phosphorus. No- 
tice in the first place the excellent agreement of the theo- 
retical predictions with the available experimental data. 
In the limit of pure single-particle behavior, the points 
should lie on straight lines joining the neutron closures 
N=20, N=28 and N=32, corresponding to the ones in 



Fig. [2] for the potassium chain. This is clearly not the 
case, because of the coupling of the proton-single parti- 
cle states with the neutron excitations. As we can gather 
from the figure, the splitting in the potassiums drops 
rapidly when two neutrons are added to the N=20 clo- 
sure, meaning that the states have not a single-particle 
character anymore. At N=28 the single-particle behav- 
ior is recovered and the inversion of the two states that 
takes place is dictated by the effective single-particle en- 
ergies. The crossing of the two orbits just reflects the 
fact that the neutron-proton interaction is more attrac- 
tive between the orbits Ofi and Oda than between the 

2 2 

orbits Of? and Isi. When neutrons flll the Ipa orbit, 

22 2 

the opposite is true and the ground-state spins of 
and ^^K bounce back to The situation is less clear 
cut in the chlorines, because its three proton holes pro- 
duce enhanced mixing. This can be seen also in Fig. [3l 
at N=20, ■^^Cl has the normal ordering of spins that cor- 
responds to a ground-state proton configuration (Isi)^ 

(Od|)^ J= 1^ with an excited i"*" state corresponding to 
the configuration (Isi)^ (Oda)^. Increasing the number 
of neutrons reduces the Oda-lsj. splitting and increases 
the mixing among different proton configurations. Al- 
ready at N=:26 both states approach the configuration 
(Is^)^ (Od|)^ and become nearly degenerate as the pro- 
tons enter in a pseudo-SU3 regime. The heavy chlorines 
may be viewed as one proton coupled to the correspond- 
ing (A-1) sulfur isotope. At N=26, "^^S is deformed and 
at N=28 the neutron closure has vanished in ^^S; this ex- 
plains why, in the chlorine chain, there is no trace of the 
single-particle effects that are so visible in the potassium 
chain when we cross N=28. Mutatis mutandis, these ar- 
guments hold for the phosphorus isotopes. 



VII. CONCLUSIONS 

In summary, we have produced a new effective interac- 
tion, SDPF-U, suitable for unrestricted Ohio large-scale 
shell- model calculations in the sd — pf valence space, for 
nuclei with proton number ranging from Z=8 to Z=20 
and neutron number from N=20 to N=40. We have pre- 
sented a first series of theoretical predictions and have 
compared them with the experimental results. Special 
attention has been devoted to the behavior of the neutron 
rich silicon and magnesium isotopes and to the vanishing 
of the N=28 neutron shell closure in ''^Si and '"'Mg. We 
have also explored the evolution of the energy splitting 
between the lowest ^ and | states in the neutron-rich 
isotopes of potassium, chlorine and phosphorus, finding 
an excellent agreement with the experimental results. 
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TABLE V: Centroids of the SDPF-U(*) and SDPF-NR in- 
teractions for A=18 (in MeV). The standard mass depen- 
dence is adopted in the calculations. The sin- 
gle particle energies (in MeV) on a core of ^^O are the fol- 



lowing 


: 0d5/2=-3.70; lsi/2=- 


2.92; 0d3/2 = l-90; 


017/2=6.22; 


lP3/2 = 


=6.31; 0f5/2 = 


=11.45; lpi/2 


=6.48. 




jl 




T 


SDPF-U 


SDPF-NR 


0d5/2 


Od5/2 





-3.19 


-3.19 


0d5/2 


0d5/2 


1 


-0.53 


-0.61 


Ud5/2 


lSi/2 





-o.iz 




0d5/2 


ISi /o 


1 


0.14 


0.14 


0d5/2 


Od3 /9 





-3.74 


-3.74 


0d5/2 


Odo/o 

" o/ z 


1 


-0.33 


-0.29 


ISl /2 


lSx/2 





-3.51 


-3.51 


lSl/2 




1 


-2.15 


-2.15 


ISl /2 


0d3/2 





-3.01 


-3.01 


ISl /2 




1 


0.00 


0.00 


0d3/2 


0*^3/2 





-2.67 


-2.67 


0d3/2 


1 2, 


1 


-0.31 


-0.40 


0d5/2 


0f7/2 





-2.55 


-2.85 


Od5/2 


0f7/2 


1 


0.00 


0.00 


0d5/2 


lP3/2 





-1.51 


-2.50 


0d5/2 


lP3/2 


1 


-0.31 


-0.08 


0d5/2 


0f5/2 





-3.17 


-4.16 


0d5/2 


0f5/2 


1 


-0.32 


-0.09 


0d5/2 


lPl/2 





-1.77 


-3.25 


0d5/2 


lPl/2 


1 


-0.48 


-0.08 


lSl/2 


0f7/2 





-1.89 


-1.89 


lSl/2 


0f7/2 


1 


-0.30 


-0.30 


lSl/2 


lP3/2 





-1.44 


-1.07 


lSl/2 


lP3/2 


1 


-1.12 


-1.41 


lSl/2 


0f5/2 





-2.48 


-3.17 


lSi/2 


0f5/2 


1 


-0.05 


0.18 


J-Sl/2 


1Di /o 





-i.yi 


-4.zi 


lSl/2 


lPl/2 


1 


0.05 


0.43 


0d3/2 


0f7/2 





-2.43 


-2.43 


0d3/2 


0f7/2 


1 


-0.88 


-0.88 


0d3/2 







-1.72 


-2.26 


0d3/2 


lP3/2 


1 


-0.30 


-0.15 


0d3/2 







-2.39 


-3.08 


0d3/2 




1 


0.27 


0.50 


0d3/2 


1Di /I 





-2.05 


-2.79 


0d3/2 


lpl/2 


1 


0.03 


0.41 


Of7/2 


0f7/2 


1 


-0.27 


-0.37 


Of7/2 


lP3/2 


1 


0.15 


0.12 


0f7/2 


0f5/2 


1 


0.03 


-0.04 


0f7/2 


lPl/2 


1 


0.20 


0.07 


lP3/2 


lP3/2 


1 


-0.10 


-0.81 


lP3/2 


0f5/2 


1 


0.36 


-0.16 


lP3/2 


lPl/2 


1 


0.16 


-0.62 


0f5/2 


0f5/2 


1 


0.12 


-0.07 


Of5/2 


lPl/2 


1 


0.14 


-0.05 


lPl/2 


lPl/2 


1 


0.25 


-0.46 



(*) without the global monopole corrections (see text). 
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